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Abstract 

The theory of strong interactions, QCD, is described in terms of a few parameters, namely the strong 
coupling constant Us and the quark masses. We show how these parameters can be determined 
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scaling method, which allows to trace the energy dependence of Os and quark masses over several 
orders of magnitude. We also discuss methods designed to reduce the effects of finite lattice spacing 
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The theory of strong interactions, QCD, is described in terms of a few parameters, namely 
the strong coupling constant as and the quark masses. We show how these paraineters can 
be determined reliably using computer simulations of QCD on a space-time lattice, and by 
employing a finite-size scaling method, which allows to trace the energy dependence of and 
quark masses over several orders of magnitude. We also discuss methods designed to reduce 
the effects of finite lattice spacing and address the issue of computer resources required. 

1 The Standard Model of particle physics 

Over the last few decades, particle physicists have explored the fundamental forces down 
to distance scales of « 10~^* m. It was found that the experimental observations are de- 
scribed to very high accuracy by a theory which is known as the Standard Model of particle 
physics. During the 1990s in particular, the predictions of this theoretical framework have 
been put to very stringent tests in accelerator experiments across the world. Perhaps one 
of the most impressive examples of its predictive power, the line shape of the Z-resonance 
in e^e~ scattering, is shown in Fig. [l|. 

The Standard Model describes the interactions of the fundamental constituents of mat- 
ter through electromagnetic, weak and strong forces in terms of three different quantum 
gauge theories. The success of the Standard Model is not only a consequence of the math- 
ematical simplicity of its basic equations, but also because the forces they describe are rela- 
tively weak at the typical energy transfers in current experiments of about 10 — 100 GeV. |^ 
The strengths of the interactions are characterized by so-called coupling constants. When 
the forces are weak, the predictions of the theory can be worked out in terms of an expan- 
sion in powers of these coupling constants, a procedure known as perturbation theory. For 
instance, in Quantum Electrodynamics (QED), the quantum gauge theory describing the 
interactions between electrons and photons, the coupling constant is the well-known fine 
structure constant a ~ 1/137. Its smallness guarantees that only a few terms in the power 
series are sufficient in order to predict physical quantities with high precision. 

The gauge theory for the strong force is called Quantum Chromodynamics (QCD), in 
which quarks and gluons assume the roles of the electrons and photons of QED. Quarks 
are the constituents of the more familiar protons and neutrons. The coupling constant of 
QCD, as, then characterizes the strength of the interaction between quarks and gluons in 

"In particle physics it is customary to use "natural units" where the speed of light, c and Planck's constant, h 
are set to one and energies as well as masses are given in GeV. As an orientation note that rriproton ~ 1 GeV, 
where 1 GeV = 1.602 ■ 10"'^ J. 
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a similar manner as the fine structure constant does in QED. One important property of all 
coupling "constants" in the Standard Model is that they depend on the energy transfer jj, 
in the interaction process. In this sense they are not really constant, and one usually refers 
to them as couplings that "run" with the energy scale. For instance, at ^ « 100 GeV the 
strong coupling constant has been determined as as = 0.12. Although this is much larger 
than the fine structure constant of QED, the method of perturbation theory still works well. 
However, if the energy scale /i is decreased from 100 GeV it is found that the value of 
as increases. In fact, at yu « 1 GeV it becomes so large that perturbation theory cannot 
be relied upon any more. It is then obvious that particle theorists require a tool which is 
able to deal with large values of ag. In other words, what is needed is a non-perturbative 
method to work out the predictions of QCD in this situation. 

As mentioned before, the simple and elegant theory of QCD is formulated in terms 
of quarks and gluons. Evidence for their existence has been accumulated in scattering 
experiments at high energies. Yet what is observed in experiments at low energies, say, fj, < 
1 GeV are protons, neutrons, 7r-mesons and many other particles, all known as hadrons. 
In fact, a striking property of QCD is "confinement", which means that quarks and gluons 
cannot be produced in experiments. Intuitively this property of QCD can be understood in 
terms of the strong growth of forces between quarks as their separation increases. Thus, 
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the only directly observable particles are the bound states of quarks and gluons, i.e. the 
hadrons. 

Particle physicists are then faced with the task of connecting the theoretically rather 
simple regime of QCD at high energies with the properties of protons, 7r-mesons and other 
hadrons observed at low energies. This task is made all the more difficult since analytic 
methods such as perturbation theory fail completely to describe the world of hadrons. 
This has led to the development of numerical techniques, namely computer simulations 
of QCD formulated on a discrete lattice of space-time points. This method allows for a 
non-perturbative treatment of the theory in the low-energy regime. 

The basic reason why the two regimes can be connected is the fact that QCD contains 
only the fundamental gauge coupling as and the masses of the quarks as free parameters. 
All observables, such as the mass of the proton, can in principle be predicted in terms of 
these quantities. The ALPHA Collaboration has embarked on a project to connect low and 
high energies in practice, byjneans of extensive computer simulations of lattice QCD. Fol- 
lowing the original idea of B, new methods have been developed and applied, which have 
yielded rather precise results. It is then possible to turn the tables: starting from the accu- 
rate experimental information on the properties of hadrons, the fundamental parameters of 
QCD can be determined numerically. 

In particular, the high-energy behaviour of the strong coupling constant as is given by 



a^ oc 



(1) 



This implies that the energy dependence of as is specified completely in terms of a single 
parameter. A, measured in energy units. It is only natural to take A as a basic parameter of 
QCD instead of the energy-dependent as- The ratio of A to the proton mass is computable 
in lattice simulations - but not in perturbation theory. Thus, one of the main goals of the 
ALPHA Collaboration is the precise determination of A. 

With QCD being one of the pillars of the 
Standard Model, it is clear that the precise 
knowledge of its parameters such as A is impor- 
tant for the ongoing quest for generalizations of 
the Standard Model. Such more complete theo- 
ries are needed to describe the early stages of the 
universe and are also expected to be relevant at 
energies which will be accessible at future par- 
ticle colliders. 



L 



2 Lattice QCD 
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Figure 2. Two dimensional slice of a 4- 
dimensional space-time lattice. Red (yellow) 
points at the bottom (left) are identified with 
those at the top (right) by periodic boundary 
conditions. 



In the mathematical formulation of QCD, the 
basic quantities are quark and gluon fields, 
which are functions of the space-time coordi- 
nates Xfj,, with xq identified with time. The clas- 
sical field equations, which describe their dy- 
namics, are differential equations - generalizations of Maxwell's equations for electromag- 
netism. To allow for a numerical treatment, it is then natural to discretize the differential 
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equations with a discretization length a, termed the lattice spacing. This procedure turns 
differential operators into finite difference operators and the fields are defined only at the 
points of a space-time lattice, illustrated in Fig. |[ 

Quantization is achieved by Feynman's path integral representation. It involves inte- 
grations over all degrees of freedom weighted with the exponential of the classical action. 
Let denote an observable, represented e.g. by a combination of quark and antiquark 
fields. Its expectation value, {ft), is defined as 

(17) ^ 1 J I?[[/]I?[v;,V]f^e-^'5[^l~^^[^''^''''l, (2) 

where Z is fixed by the condition (1) = 1. To further prepare for an evaluation of the path 
integral on a computer, the quark degrees of freedom are integrated out analytically. The 
expression for (il) then becomes 

(f^) = i y V[U] {dct(I?[t/])}^'' e-^^^^l (3) 

fJeff denotes the representation of ft in the effective theory, where only gluon fields remain 
in the path integral measure. Equation (|^) requires some further explanation: 

• D[U] denotes the Dirac operator, and for simplicity of presentation we have consid- 
ered QCD with Ni quarks of equal mass rn, which enters D[U]. For hadron physics 
at small energies the case A't = 3 is most relevant. The heavier quarks decouple from 
the dynamics up to small effects; 

• The gluon field is represented by link variables U{x, n) connecting the sites x and x + 
ajl (shown as _^ in Fig.^). The measure V[U] is the product of integration measures 
for each link; 

• The action S'g[C^] is a sum of local terms over the lattice, coupling only gluon vari- 
ables within one plaquette ( in Fig. By contrast, the effective interaction re- 
sulting from the integration over the anticommuting quark fields is of infinite range: 
although the finite difference Dirac-operator D[J7] is local, <let{D\U]) couples gluons 
at arbitrary distances. 

• Using the representation eq. (^, our expression for {Vl) has the form of a thermal 
average in statistical mechanics, and one may use stochastic sampling to evaluate the 
integral if the space-time volume is made finite by restricting < a;^; < L. 

• As an example, consider ($(a;)$t (y)) where $(a;) is a suitable combination of quark 
fields at a point x, which has the quantum numbers of some hadron. The expectation 
value is then proportional to the quantum-mechanical amplitude for the propagation 
of the hadron from point y to x, from which the mass of the hadronic bound state may 
be obtained. 

The lattice formulation sketched above leads to a mathematically well-defined expression 
for (il). In particular, the typical infinities which are encountered in field theoretical ex- 
pectation values are absent. 

The exact treatment of the determinant, eq. (||), still presents a major challenge, even 
on today's massively parallel computers. In many applications one has therefore set A^f = 
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0. This - very drastic - approximation defines the so-called quenched approximation. 
Physically it means that the quantum fluctuations of quarks are neglected and only those 
due to the gluons are taken (exactly) into account. Although it turns out that the quenched 
approximation works well at the 10%-level, the proper treatment for A^f > is perhaps the 
most important issue in current simulations. 

Whether or not the quenched approximation is employed, one always has to address 
the problem of lattice artefacts. Let R denote a dimensionless observable, such as a ratio 
of hadron masses. Then its expectation value on the lattice differs from the value in the 
continuum by corrections of order a^: 

^lat ^ ^cont + 0(aP)^ (4) 

where the power p depends on the chosen discretization of the QCD action. The correction 
term for typical values of a can be quite large, and an extrapolation to the continuum limit 
is then required to obtain the desired result. If p is large, the rate of convergence to the 
continuum limit is high, and hence the extrapolation is much better controlled. 

Let us now return to the problem of determining the energy dependence of quantities 
such as the running coupling a^ . Obviously the numerical simulation is possible only if the 
number of points of the lattice, (L/a)^, is not too large; typical lattice sizes are L/a < 32. 
Therefore, besides a also the effect of the finite value of L has to be considered. It is known 
rather well that L = 2 fm is sufficient for most quantities, and that much smaller volumes 
would lead to unacceptably large corrections to the desired L ^ oo limit. A physical box 
size of L = 2 fm together with L/a < 32 thus implies 

a > 0.05 fm. (5) 

The existence of such a lower bound on a, imposed by practicability considerations, also 
means that the energies, fx, that can be expected to be treated correctly have to satisfy 

M<a"^<4GcV. (6) 

In other words, lattice QCD is well suited for the computation of low energy properties of 
hadrons, while high energies appear impossible to reach. The ALPHA Collaboration has 
developed and applied an approach to circumvent this problem, which is described below. 

3 Running coupling and quark masses 

We are now going to outline the strategy which allows to connect the low- and high-energy 
regimes of QCD in a controlled manner. Results for the A parameter and the quark masses 
will then be presented. From now on we will drop the subscript "s" on the strong coupling 
constant a^. 

3.1 Gedankenexperiment 

From the above it is evident that the restriction of numerical lattice QCD to low energies 
is necessary to avoid finite-size effects (FSE) when working with a manageable number 
of lattice sites. For the computation of an effective coupling, this problem can be circum- 
vented since one has great freedom in the definition of such a coupling. Even the strength 
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of finite-size effects serve as a measure of the interactions of the theory, and thus a suitably 
chosen FSE may be used to define El an 



effective coupling g{L) with 



An 



fi=l/L. 



(7) 



It depends on ("runs with") the energy scale 1^1= 1/L. Such a coupling can be computed 



in the regime ^ ^ 



requiring only a moderate resolution of the space-time world; 



L/a ~ 0(10) points per coordinate are sufficient. 

In a numerical calculation, the success of this general idea will depend on a few prop- 
erties of the coupling. It must be computable with good statistical precision in a Monte 
Carlo (MC) simulation, and with small discretization errors. Furthermore, one would like 
to know its scale dependence analytically for large energies. This is achieved by determin- 
ing the perturbative expansion of the so-called /3-function, 



m) = -L^ -f {bo + hf + + ...}, 



(8) 



constant 



to the order indicated above, or even higher. 

In QCD, a coupling with 
properties could indeed be found 
Its definition starts from QCD in a box 
of size L**. Periodic boundary con- 
ditions are imposed for the fields as 
functions of the three spatial coordi- 
nates, and Dirichlet boundary condi- 
tions are set in the time-direction, as 
shown in Fig. ||. The Dirichlet bound- 
ary conditions are homogeneous except 
for the spatial components of the gluon 
gauge potentials, Ak- With this choice 
of boundary conditions the topology is 
that of a 4-dimensional cylinder 

Classically, these boundary condi- 
tions lead to a homogeneous (x and 
t independent) colour-electric field in- 
side the cylinder The walls at time 

< = and t = r act quite similarly to the plates of an electric condensor The strength 
of the QCD interactions is conveniently defined in terms of the colour-electric field at the 
condensor plates: 



classically: 
electric field 

E{x) 

all quantum fields 4>: 

(j){x + Lk) = (l){x) 
periodic in space 



II 



time t 

Figure 3. The boundary conditions used for the definition 
of the finite-size coupling g{L) in QCD. 



f{L) = 



lassical 



{E) 



(9) 



Here i? is a special colour component of the electric field. 

For weak coupling, i.e. for small L, the path integral is dominated by field configu- 
rations which correspond to small fluctuations about the classical solution. On the other 
hand, for L « 1 fm they may deviate significantly from the classical solution, and this can 
be realized in a MC-simulation of the path integral. 
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3.2 The running coupling 



The energy dependence of the coupUng g'^{L) can be determined recursively through a 
number of compute-effective steps. Figure ^illustrates the implementation of the recursion 



-g\2L)^aifiL)). 



(10) 



The function a describes the 
change in the coupling when 
the physical box size is dou- 
bled. Since the relevant en- 
ergy scale for the running 
of the coupling (/i = l/L) 
is separated from the lattice 
spacing a, one may com- 
pute a recursively over sev- 
eral orders of magnitude in 
fi, whilst keeping the num- 
ber of lattice sites at a man- 
ageable level, i.e. L/a = 
0(10). 

The results of one 
horizontal step in Fig. ^ 
performed on the APE- 
computers at NIC/Zeuthen, 
is displayed in Fig. ^ They 
are shown as a function of 
the resolution a/L. Because 
the coupling and the details 
of the discretization were 

chosen with great care, the dependence on the resolution is tiny and the lattice numbers 
can be extrapolated to the continuum a/L ^ 0. To arrive at the continuum limit, each of 
the horizontal steps in the figure is indeed repeated several times with different resolutions ! 

In this way, the continuum func- 
tion (j{u) was obtained for a range 
of M = g^{L). Starting from a min- 
imal value of u, the recursive appli- 
cation of eq. (|l^ yields the points 
in Fig. ^. They may be compared 
to the integration of the differential 
equation towards low values of 
fi, starting at the lowest value of a. 
One then makes the important ob- 
servation, that eq. truncated at 
the indicated order, is quantitatively 

verified in the region of low enough a. It may hence also be used to compute the A- 



Figure 4. The recursion to compute g{L) with L = Lq ^ 2Lq — ► 
ALq — 8-Lo- It is iterated further until one arrives at 2^Lq, in practice 
n ~ 8. One always alternates between keeping a fixed and increasing 
L 2L and keeping L fixed and increasing a. 
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Figure 5. The extrapolation of the function rp[2L) = 
(j{g'^{L)) from the data at finite lattice spacing a to the con- 
tinuum (blue square). 
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Figure 6. The energy dependence of the coupHng as computed from lattice QCD (points) compared to the per- 
turbative expression eq. (H) with the f)2-term (full line) and without. 



parameter from (g = g{L = l//i)) 
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1 1 bi 



(11) 



with negligible errors due to higher-order terms that are not included. 

The attentive reader will have noticed that the computation explained so far, "only" 
determines the dependence of a on the combination /i/A, while in the figure we show it 
as a function of // in physical units for A^f = 0, i.e. in the quenched approximation. The 
missing link is to connect the lowest energy = 1/ L,„ax contained in the figure, to a low 
energy - experimentaJly accessible - property of a hadron. The connection with the decay 
rate of the K-meson u, which we do not describe here, yields the physical units shown in 
the figure as well as the final result 

A = 238(19)McVfor^f = 0. (12) 

Currently, the ALPHA Collaboration is extending these computations to the numerically 
very demanding case of TVf = 2. First promising results, which are still awaiting a final 
check for the absence of discretization errors are shown on the rh.s. of Fig. ^ The step to 
connect imax to a low-energy observable has yet to be performed, and hence the results are 
plotted as a function of /i/A. It is worth pointing out that, in order to ensure efficiency and 
correctness of thesai^ > simulations, the development and testing of MC-algorithms is 
very important Eliillij. 



3.3 Quark masses 

The masses of the quarks are of a very different nature from the mass of the electron. The 
fundamental difference is the property of confinement mentioned in section 1 . It means 
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m(2 GeV)[MeV] 



light strange charm beauty 



10 100 1000 
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Figure 7. Quark masses determined by the ALPHA collaboration in the A^f = approximation E3tjl 
The label "light" denotes the average mass pfithe up and the down quark mass. Green shaded areas show the 
error band quoted by the particle data group ESI- 



that a quark can not be prepared in isolation to perform an experimental measurement of 
its mass. Consequently, quark masses have to be understood as parameters of the theory, 
and their proper definition very much resembles that of an effective coupling. In particular, 
they have an energy dependence similar to the one of a. When defined in a natural way E3, 
the energy dependence is the same for all quark flavours; in other words, ratios of quark 
masses do not depend on the energy. The overall energy dependence has been determined 
with the jSame strategy as the one described for a and in fact with similar precision (see 
Fig. 2 in EJ). 

The numerical values of the quark masses (at a particular low energy scale) are con- 
veniently extracted by relating them to the masses of quark-antiquark bound states with 
pseudoscalar quantum numbers: tt-, K-, D-, B-mesons. Results of the ALPHA Collabora- 
tion obtained in the A'f = approximation are shown in Fig. ^in two different (common) 
conventions. ^ The upper row shows the running masses m in the so-called MS scheme at 
a renormalization energy — 2 GeV, while the lower one shows the so-called renormal- 
ization group invariant quark masses. The latter are related to the running quark masses 
via 



M = lim m{[i) 2bo {g{n)) 



-do/2fco 



(13) 



where do characterizes the asymptotic behavior of m(/x) for large energy. It is computable 
in perturbation theory and has a value of do = 8/ (47r)^. 



4 Improvement 

Equation (Q) shows that lattice and continuum quantities differ by discretization errors. In 
Wilson's formulation of lattice QCD these discretization effects are of order a, which can 
be quite large. People have therefore tried to find an improved version of the Wilson action, 
so that the deviation of observables computed on the lattice to their continuum values is of 
O(a^), or even higher. 

''The ratio of strange to light quark masses actually relies on 0. In the special case of the b-quark, whose mass 
is large compared to the QCD scale, A, an expansion in terms of 1/ has been used and only the lowest order 
term was considered tj. 



9 



In order to illustrate the problem, let us consider the solution of differential equations 
by means of numerical integration. Given the initial value problem 

:7- = /(a;,2;), y{xo)^yo, (14) 

ax 

the well-known Euler method provides a formula for its numerical solution, i.e. 

yj+i = % +«/(^j>%)> j = 0,1,2,.... (15) 

Here the quantity a denotes the finite step size used to approximate the derivative dy/dx in 
terms of finite differences. In order to compute the solution at point x = xo+na for a fixed 
value of a one has to perform n iteration steps. The central question then is by how much 
this solution deviates from the exact one. For the simple Euler method one can show that 
this so-called truncation error is proportional to a^. Obviously the accuracy of the solution 
will be larger for small a, but then the number of iterations, n, may become very large. 
It is much more efficient to use an improved solution scheme, such as the Runge-Kutta 
method. Here the function / in eq. ( |l5| ) is replaced by a more complicated expression, 
which is chosen in such a way that the truncation error is proportional to a^. 

The analogy of this example with lattice QCD is obvious: the step size is now the 
lattice spacing, and the lattice action must be chosen such that the leading discretization 
error is of higher order in a. However, in a quantum field theory like QCD the search for 
an improved action is not straightforward: although one can quite easily write down an 
action which formally has O(a^) artefacts, the interactions between quarks and gluons in 
the discretized theory can again introduce 0(a) errors at the quantum level. 

To be moie explicit, let 5w denote the Wilson action for lattice QCD. Sheikholeslami 
and Wohlert 113 have shown that it is sufficient to add just one more interaction term, cor- 
responding to an action 6S, to 5w in order to cancel its leading discretization error pro- 
portional to a. To ensure that quantum effects do not introduce 0(a) effects "via the back 
door", the term SS has to be multiplied by a coefficient Cgw This improvement coefficient 
must be suitably tuned in order to achieve the complete cancellation of 0(a) artefacts at 
the quantum level for observables such as hadron masses. The 0(a) improved version of 
Wilson's action then reads 

•S'sw = Sw + CswSS. (16) 

The ALPHA CoUabyotatitrnJias developed and applied a method to determine Cgw in com- 
puter simulations E2lEill23'Ej. The method amounts to computing the expectation value of 
a pure lattice artefact as a function of Cgw The value of Csw at the point where the artefact 
vanishes then defines the improved action for the corresponding value of the lattice spac- 
ing. Thus, non-perturbative methods are not only used to compute physical observables, 
but also to improve the reliability of the numerical treatment of the theory as such. The 
improved action is only slightly more complicated to implement in practical simulations. 
The increased effort is easily offset by the fasteu. convergence to the continuum limit. As a 
result, the action with as determined in E3l2j has by now become a standard for preci- 
sion lattice QCD computations. In particular this improvement was essential in obtaining 
most of the results quoted in the previous section. 

Of course, one also has to verify that observables computed with the improved action 
indeed approach the continuum limit with a tale proportional to a^. Tests of this kind have 
been performed successfully (see, e.g. ref. lJ), and other examples for the effectiveness 
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of the improvement programme for the evaluation of many hadronic observables can be 
found in u. 



5 Machines and resources 

Although conceptual advances like 0(a) improvement are of great importance, progress in 
lattice QCD is also dependent on the availability of sufficient computer resources. 

In the 1980s most simulations were performed on vector supercomputers like the Cray 
X-MP. Since then the demand for increased performance has led to the development of 
massively parallel machines, which are now widely used. Lattice QCD is a problem which 
lends itself easily to parallelization; the total volume can be divided into many sublattices, 
which are distributed over a grid of processors. The latter can perform the same task on 
independent data. Furthermore, the communication pattern for lattice QCD is simple, since 
most algorithms only require nearest-neighbour communications. 

The AT PHA CnUahnraunn has mostly used parallel computers from the APE family 
of machines Eju3€3l3. The latest generation, APEmille, has been developed jointly by 
INFN and DESY. The smallest entity in the APEmille processor grid consists of a cube of 
2x2x2 processors. These cubes are then connected to form larger grids. Despite its 
very conservative clock speed of 66 MHz, each processor achieves a peak performance of 
528 MFlops, ^ thanks to the optimization for complex arithmetic: the operation a x b + 
c, which requires 8 floating point operations for the three 32-bit complex numbers a, b 
and c, is performed in one clock cycle. The total peak speed of the current installation of 
APEmille machines at DES Y-Zeuthen amounts to more than 500 GFlops (single precision). 

The programming language for APEmille is TAO, which has a FORTRAN-like syntax, 
but also includes special features designed to facilitate parallelization and coding, and al- 
lows to achieve a high proportion of the peak speed. One such feature is the easy access 
to any desired number of the 512(!) registers. The typical efficiency of ALPHA-programs 
on APEmille is about 30% of the peak speed. A similar figure has been achieved by 
other lattice QCD collaborations on machines like the-.Cray T3E, but at the price of hav- 
ing to code the core routines in assembler language Ej. An interesting figure-of-merit is 
the price/performance ratio. For APEmille it is (8 — 10) Euro per sustained MFlops and 
there are efforts under way to reduce this number much further In addition it was also 
demonstrated that PC clusters with Myrinet network can achieve about 3 — 4 Euro/MFlops 
sustained E^l. These developments let us await the future with optimism. 



6 Future 

Moving towards the realistic case of iVf = 3 should soon be possible. Once this has been 
achieved, the most precise determinations of the fundamental parameters of QCD may 
come from the low energy hadron spectrum combined with lattice QCD to evaluate the 
theory. The methods developed in this project are expected to play an important role in 
this program. 

Furthermore these methods and related ones Ej>E31±3 will improve the reliability of the 
determination of properly normalized weak decay (and mixing) amplitudes of hadrons c3. 

'^The unit 1 MFlops denotes one million floating point operations per second. 
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Their knowledge is again of vital importance in tests of the Standard model and the search 
for an even more fundamental theory. 



The results presented in this article were obtained on the APE 100 and APEmille instal- 
lations at Zeuthen. The total CPU time required was about lO*^ processor-hours on APEIOO 
and 0.2 • 10^ processor-hours on APEmille, respectively. 
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